INCOMPRESSIBLE LIMIT OF THE COMPRESSIBLE 
NON-ISENTROPIC MAGNETOHYDRODYNAMIC EQUATIONS 
WITH ZERO MAGNETIC DIFFUSIVITY 
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Abstract. We study the incompressible limit of the compressible non- 
isentropic magnetohydrodynamic equations with zero magnetic diffusivity and 
general initial data in the whole space R** (d = 2,3). We first establish the 
existence of classic solutions on a time interval independent of the Mach num- 
ber. Then, by deriving uniform a priori estimates, we obtain the convergence 
of the solution to that of the incompressible magnetohydrodynamic equations 
as the Mach number tends to zero. 



1. Introduction 

This paper is concerned with the incompressible Umit to the compressible non- 
iscntropic magnetohydrodynamic (MHD) equations with zero magnetic diffusivity 
and general initial data in the whole space Mf^ {d ~ 2, 3). 

In the study of a highly conducting fluid, for example, the magnetic fusion, it 
is rational to ignore the magnetic diffusion term in the MHD equations since the 
magnetic diffusion coefficient (resistivity coefficient) is inversely proportional to the 
electrical conductivity coefficient, sec [10]. In this situation, the system, describing 
the motion of the fluid in M'', can be described by the following compressible non- 
isentropic MHD equations with zero magnetic diffusivity: 

+ div(pu) = 0, (1.1) 

dt{pn) + div (pu (g) u) + Vp = (V X H) X H + div^-, (1.2) 

^tH- V X (u X H) = 0, divH = 0, (1.3) 

dt£ + div {vi{£' + p)) = div((u X H) X H) + div(u1' + nVO). (1.4) 

Here p denotes the density, u e the velocity, H e M'' the magnetic field, and 6 
the temperature, respectively; is the viscous stress tensor given by 

^ = 2^D(u) + Adivu I<j 

with D(u) = (Vu + Vu^)/2, and 1^ being the c? x d identity matrix, and Vu^ the 
transpose of the matrix Vu; £ is the total energy given hy £ = £' + |Hp/2 and 
£' = p (e + |up/2) with e being the internal energy, p|up/2 the kinetic energy, and 
|Hp/2 the magnetic energy. The viscosity coefficients A and ji of the fluid satisfy 
2/^ + dA > and > 0; k > is the heat conductivity. For simplicity, we assume 
that p, A and k, are constants. The equations of state p — p{p, 8) and e = e(p, 9) 
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relate the pressure p and the internal energy e to the density p and the temperature 
9 of the flow. 

For the smooth solution to the system (1.1)-(1.4), we can rewrite the total energy 
equation (f .4) in the form of the internal energy. In fact, multiplying (1.2) by u 
and (1.3) by H, and summing the resulting equations together, we obtain 



- (^-p|u|^ + -|H|^ j + -div(p|u|^u) + Vp • u 

= div* • u + (V X H) X H • u + V X (u X H) • H. (1.5) 

Using the identities 

div(H X (V X H)) = |V X Hp - V X (V X H) • H, (1.6) 

div((u X H) X H) = (V X H) X H • u + V X (u X H) • H, (1.7) 
and subtracting (1.5) from (1.4), we thus obtain the internal energy equation 

dt{pe) + div(/ciue) + (divu)p = * : Vu + kM, (1.8) 
where ^' : Vu denotes the scalar product of two matrices: 

* ^ = E § (1^ + IJ)' + = 2m|B(u)|2 + A(trD(u))2. 



Using the Gibbs relation 



edS ^de+pd\-] , (1.9) 



we can further replace the equation (1.8) by 

dt{pS)+dw{pSn)^'^ -.Vu + kM, (1.10) 

where S denotes the entropy. 

In the present paper, we assume that k = in (1.10). Now, as in [27], we 
reconsider the equations of state as functions of S and p, i.e., p — R{S,p) and 
6 = Q{S,p) for some positive smooth functions R and Q defined for all S and 
p > 0, and satisfying dR/dp > 0. For instance, we have p = p^l'^e^^l'^ for ideal 
fluids. Then, by utilizing (1.1) together with the constraint divH = 0, the system 
(1.1), (1.2), (1.4) and (1.10) can be rewritten as 

A(S,p)(dtP + (u ■ V)p) + divu = 0, (1.11) 

R{S,p)(dt\x + (u • V)u) + Vp = (V X H) X H + div*, (1.12) 

- curl (u X H) = 0, divH = 0, (1.13) 

R{S,p)Q{S,p){dtS + (u ■ V)5) = : Vu, (1.14) 

where 

Considering the physical explanation of the incompressible limit, we introduce 
the dimensionless parameter e, the Mach number, and make the following changes 
of variables: 

pix, t) = p'(a;, et), S{x, t) = S'{x, ei), 
u(a;, t) = eu'^{x, et), H(a;, t) ~ eH'(x, et), 
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and 

As the analysis in [27], we use the transformation p'^{x,et) = pe^'' for some 
positive constant p. Under these changes of variables, the system (1.11)-(1.14) 
becomes 

a'{S', eq'){dtq' + {u' ■ V)q'') + -divu' = 0, (1.16) 

e 

r'{S',€q'){dtu'' + (u' • V)u') + -Vq" = (curlH^) x + div^'^ (1.17) 

e 

dtH" - curl (u' X H"^) = 0, divH' = 0, (1.18) 
b''{S\eq'){dtS'' + {u' ■ \7)S') = e^^' : Vu% (1.19) 
where we have used the abbreviations 'I''^ = 2/i'^D(u'^) + A'^divu'^ 1^ and 

a%S\ eq^ := A{S\ pe^'^' )pe^'^' = „,g'^\^.. • ^ , (1.20) 

- - R{S'^,pe'^'' ) OS s=pe"i' 

r\S\ eq^ := 'g. \ b^iS\ eq^) := R{S^eq^)QiS\eq^). (1.21) 

Formally, we obtain from (1.16) and (1.17) that Vg^ — ?► and divu' = as 
e — > 0. Applying the operator curl to (1.17), using the fact that curlV = 0, and 
letting e — !■ and yu'^ — s- /i > 0, we obtain 

curl (r(5, 0)i^t^^ + v Vv) - (curlH) x H - /iAv) = 0, 

where we have assumed that (5"^, q'^, u^, H*^) and r'^(S"^, eq"^) converge to (5*, 0, v, H) 
and r{S,0) in some sense, respectively. Finally, applying the identity 

curl (u x H) = u(divH) - H(divu) + (H • V)u - (u • V)H, (1.22) 

we expect to get the following incompressible non-isentropic MHD equations 

r(5,0)(atv + (v-V)v)-(curlH) x H + Vtt = ^iAv, (1.23) 

atH + (v-V)H-(H-V)v = 0, (1.24) 

dtS + {v-\7)S ^0, (1.25) 

divv = 0, divH = (1.26) 

for some function tt. 

The aim of this paper is to establish the above limit process rigorously in the 
whole space M.'^. 

Before stating our main results, we review the previous related works. We begin 
with the results for the Euler and Navier- Stokes equations. For well-prepared initial 
data, Schochet [33] obtained the convergence of the compressible non-isentropic 
Euler equations to the incompressible non-isentropic Euler equations in a bounded 
domain for local smooth solutions. For general initial data, Metivier and Schochet 
[27] proved rigorously the incompressible limit of the compressible non-isentropic 
Euler equations in the whole space M''. There are two key points in the article 
[27]. First, they obtained the uniform estimates in Sobolev norms for the acoustic 
component of the solutions, which are propagated by a wave equation with unknown 
variable coefficients. Second, they proved that the local energy of the acoustic 
wave decays to zero in the whole space case. This approach was extended to the 
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non-isentropic Euler equations in the exterior domain and the full Navier-Stokes 
equations in the whole space by Alazard in [1] and [2], respectively, and to the 
dispersive Navicr-Stokcs equations by Lcvermore, Sun and Trivisa [26]. For the 
spatially periodic case, Mctivicr and Schochet [28] showed the incompressible limit 
of the one-dimensional non-isentropic Euler equations with general data. Compared 
to the non-isentropic case, the treatment of the propagation of oscillations in the 
isentropic case is simpler and there are many works on this topic. For example, see 
Ukai [35], Asano [3], Desjardins and Grenier [7] in the whole space; Isozaki [16,17] on 
the exterior domain; Iguchi [15] on the half space; Schochet [32] and Gallagher [11] 
in a periodic domain; and Lions and Masmoudi [30], and Desjardins, et al. [8] in 
a bounded domain. Recently, Jiang and Ou [22] investigated the incompressible 
limit of the non-isentropic Navier-Stokes equations with zero heat conductivity and 
well-prepared initial data in three-dimensional bounded domains. The justification 
of the incompressible limit of the non-isentropic Euler or Navier-Stokes equations 
with general initial data in a bounded domain or a multi-dimensional periodic 
domain is still open. The interested reader can refer to [5] for formal computations 
on the case of viscous polytropic gases and [4, 28] for some analysis on the non- 
isentropic Euler equations in a multi-dimensional periodic domain. For more results 
on the incompressible limit of the Euler and Navier-Stokes equations, please see the 
monograph [9] and the survey articles [6,31,34]. 

For the isentropic compressible MHD equations, the justification of the low Mach 
limit was given in several aspects. In [23], Klainerman and Majda first studied the 
incompressible limit of the isentropic compressible ideal MHD equations in the 
spatially periodic case with well-prepared initial data. Recently, the incompressible 
limit of the isentropic viscous (including both viscosity and magnetic diffusivity) of 
compressible MHD equations with general data was studied in [14, 18, 19]. In [14], 
Hu and Wang obtained the convergence of weak solutions of the compressible vis- 
cous MHD equations in bounded, spatially periodic domains and the whole space, 
respectively. In [18], the authors employed the modulated energy method to verify 
the limit of weak solutions of the compressible MHD equations in the torus to the 
strong solution of the incompressible viscous or partial viscous MHD equations (the 
shear viscosity coeSicient is zero but the magnetic diffusion coefficient is a positive 
constant). In [19], the authors obtained the convergence of weak solutions of the 
viscous compressible MHD equations to the strong solution of the ideal incompress- 
ible MHD equations in the whole space by using the dispersion property of the wave 
equation if both the shear viscosity and the magnetic diffusion coefficients go to 
zero. 

For the full compressible MHD equations, the incompressible limit in the frame- 
work of the so-called variational solutions was studied in [24,25,29]. Recently, the 
authors [20] justified rigourously the low Mach number limit of classical solutions 
to the ideal or full compressible non-isentropic MHD equations with small entropy 
or temperature variations. When the heat conductivity and large temperature 
variations are present, the low Mach number limit for the full compressible non- 
isentropic MHD equations justified in [21]. We emphasize here that the arguments 
in [21] are different from the present paper and depend essentially on positivity of 
fluid viscosity and magnetic diffusivity coefficients. 

As aforementioned, in this paper wc want to establish rigorously the limit as 
e — > to the system (1.16)-(1.19) for /le — > /i > 0. In this case, the magnetic 
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equation is purely hyperbolic due to the lack of magnetic diffusivity. The first- 
order derivatives of H'^ in the momentum equation and magnetic equation cannot 
be controlled. It is very hard to study the system (1.16)-(1.19). To the author's 
knowledge, there is very few mathematical analysis on the system (1.16)-(1.19) 
with fixed or unfixed e, even for the iscntropic case. Our main idea is trying to 
make full use of the fluid viscosities to control the higher order derivatives of the 
magnetic field. 

Now, we supplement the system (1.16)-(1.19) with initial conditions 

(5^g^u^H^)l,=o = (5^,<?^,u^,h^) (1.27) 

and state the main results as follows. 

Theorem 1.1. Let s > d/2 + 2 be an integer. Assume that yu'^ — s- ^ > and A*^ — > A 
as e — > 0. Suppose that the initial data (5q, gg, Uq, Hq) satisfy 

\\{S'o.<1o,<^o)\\hhR') < Mo- (1.28) 

Then there exists a T > such that for any e £ (0, 1], the Cauchy problem (1.16)- 
(1.19), (1.27) has a unique solution (5%g%u%H^) e C°([0, T], iJ^(M'')), and there 
exists a positive constant N , depending only on T and Mq, such that 

||(5%g%u%ff)(i)||/f=(K.) < A^, V<e[0,T]. (1.29) 

Furthermore, if there exist positive constants S_, Nq and d such that Sq{x) satisfies 

\S'o{x)-S\<No\x\-^-', |V5^(x)| <7Vo|x|-2-^ (1.30) 

then the sequence of solutions (5*', q*^, u'^, H*^) converges weakly in L°°{0,T; H^lMf^)) 
and strongly in L'^{0,T; Hf^^(M.'^)) for all s' < s to a limit (.9, 0,v,H), where 
(5,v,H) isjhe unique solution in C([0, T], i?-'(M'^)) of (1.23)-(1.26) with initial 
data (^,v, H)|t=o (•S'oi wg, Hq), where wq G H'^(R'^) is determined by 

divwo^O, curl (r(S'o, 0)wo) = curl (r(S'o,0)vo), r(S'o, 0) lim r'(5^, 0). (1.31) 

The function tt G C([0,T] x R'^) satisfies Vtt G C([0, T], i7''-i(E^)). 

We briefly describe the strategy of the proof. The proof of Theorems 1.1 in- 
cludes two main steps: the uniform estimates of the solutions, and the convergence 
from the original scaling equations to the limiting ones. Once we have established 
the uniform estimates (1.29) of the solutions in Theorems 1.1, the convergence of 
solutions is easily proved by using the local energy decay theorem for fast waves 
in the whole space, which is shown by Metivier and Schochet in [27]. Thus, the 
main task in the present paper is to obtain the uniform estimates (1.29). For this 
purpose, we shall modify the approach developed in [27]. In fact, due to the strong 
coupling of hydrodynamic motion and magnetic field, and the lack of magnetic dif- 
fusivity, new difficulties arise in obtaining the uniform estimates for the solutions to 
(1.16)-(1.19), (1.27). First of ah, when we perform the operator {{E''}-'^L{d^))'^ 
to the continuity and momentum equations, or the operator curl to the momentum 
equations, one order more spatial derivatives arise for the magnetic field, and this 
prevents us from closing the energy estimates. Second, since the coefficients of the 
acoustic wave equations depend on the entropy, we could not get the estimates of 
l|u'^llL2(o_T;_f/=+i) directly from the system. The ideas to overcome these difficulties 
here are the following: We transfer one spatial derivative from the magnetic field 
to the velocity with the help of the special coupled way between magnetic field and 
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fluid velocity. Then, to control ||u'^||^2(g j^.^s+i^, we employ a kind of Helmholtz 
decomposition of the velocity. Third, we make full use the special structure of the 
magnetic field equation and the estimates on u to control |jH'^||£,oo(Q y./^s). 

We point out that our arguments in this paper can be modified slightly to the 
case of the the compressible non-isentropic MHD equations with infinite Reynolds 
number. We shall give a brief discussion in Section 5. 

This paper is arranged as follows. In Section 2, we give notations, recall basic 
facts, and present commutators estimates. In Section 3 we establish the uniform 
boundeness of the solutions and prove the existence part of Theorem 1.1. In Section 
4, we use the decay of the local energy to the acoustic wave equations to prove the 
convergent part of Theorem 1.1. In the last section, we consider the incompressible 
limit to the compressible non-isentropic MHD equations with infinite Reynolds 
number. 

2. Preliminary 

We give notations and recall basic facts which will be used frequently in the 
proofs. 

(1) We denote (•,■) the standard inner product in L'^{R'^) with (/, /) = 

and H'^ the usual Sobolcv space W'^''^ with norm || • in particular, || ■ ||o = || • ||. 
The notation . . . , means the summation of \\Ai\\ {i = I,-- - and 

it also applies to other norms. For a multi-index a = (ai,...,ad), we denote 
= d^l . . . d!^^ and \a\ = |q;i | -I- • • • -I- \ad\- We also omit the spatial domain R'^ 
in integrals for convenience. We use the symbols K or Cq to denote the generic 
positive constants, and C(-) and C(-) to denote the smooth functions, which may 
vary from line to line. 

(2) For a scalar function / and vector functions a, b and c, we have the following 
basic vector identities: 



div(a X b) 


= b • curl a 


- a • curlb. 


(2.1) 


V(|ap) 


= 2(a- V)a 


+ 2a X curl a. 


(2.2) 


curl (/a) 


~ f ■ curl a - 


- V/ X a. 


(2.3) 


curl (a X b) 


= (b ■ V)a - 


- (a ■ V)b + a(divb) - b(diva). 


(2.4) 


div((a X b) X c) 


~ c ■ curl (a 


X b) — (a X b) • curie. 


(2.5) 



(3) We have the following well-known nonlinear estimates [12]. 

(i) Let a = {ai,a2,ad) be a multi-index such that \a\ = k. Then, for all 
(T > 0, and f,g <E H''^"' (R.'^) , there exists a generic constant Co such that 

\\[f,d^]g\\„. <Co(||/||wi,oo||g||ff.+'.-i + ||/||ff.+.||.9||L=c,). (2.6) 

(ii) For integers k>0, l>0, k + l<a and a > d/2, the product maps 
continuously H''-''{W'^) x i7'"-'(R'') to H"-''-'- (R^) and 

\\uv\\,_k-i < K\\u\\,_k\\v\\,^i. (2.7) 

(iii) Let a > d/2 be an integer. Assume that F{u) is a smooth function such 
that F{0) = and w e H''{R'^), then F{u) e iJ'"(R'*) and its norm is bounded by 

\\F{u)\\^<C{\\uU\\u\\,, (2.8) 

where C(-) is independent of u and maps [0, co) into [0, co). 
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3. Uniform estimates 

In this section and the first part of the next section we assume that fi'^ = ^ > 
and A"^ = A for simplicity of the presentation. The general case can be treated by 
a slight modification in the arguments presented here. 

In view of [27] and the classical local existence results obtained by Vol'pert and 
Khudiaev [36] for hyperbolic-parabolic systems, the key point in the proof of the 
existence part of Theorem 1.1 is to establish the uniform estimate (1.29), which 
can be deduced from the following a priori estimate. 

Theorem 3.1. For any e > 0, let {S',q^,u',W) <= C([0, T], i7"(R'*)) be the solu- 
tion to (1.16)-(1.19). Then there exists an increasing C(-) from [0,oo) to [0,oo), 
such that 

M,{T)<Co + {T + e)C{M,{T)), (3.1) 

where 

Me{T) -.^KiTf + r llu^ll^+idr, (3.2) 
Jo 

with 

K{T):= sup ||(5^Q^u^H^)(^)|u. (3.3) 

tG[0,T] 

The remainder of this section is devoted to establishing (3.1). In the calcula- 
tions that follow, we always suppose that the assumptions in Theorem 1.1 hold. 
We consider a solution (5'%q%u%H'=) to the problem (1.16)-(1.19), (1.27) on 
C{[Q,T],H^{W^)) with initial data satisfying (1.28). 

The main idea for proving the uniform estimate (3.1) is motivated by the work 
[27] where the operator {{E'^}~'^ L{dx)y^ is introduced to control the acoustic com- 
ponents of velocity for the Euler equations. When the strong coupling of the fluid 
and magnetic filed is present, however, the arguments in [27] cannot be directly 
applied to get a uniform estimate of the acoustic parts due to lack of magnetic 
diffusion in the magnetic equation. Instead, here we transfer one order spatial 
derivative from to u*^, and then employ the fluid viscosity to control higher 
derivatives. We remark that the reason that these techniques work is due to the 
special structure of coupling between the fluid and magnetic fields. 

We begin with the estimate on the entropy 5'. 

Lemma 3.2. There exist a constant Cq > and a function C{-), independent of 
e, such that for all t G (0, T], 

WS'ml <Co + tC{M,{T)) + e^C{M,{T)). (3.4) 

Proof. For the multi- index a satisfying |a| < s — 1, denote = d"S'^. In view of 
the positivity of 6'(S"^, eg*^), we deduce from (1.19) that, 

dtfo. + iu' ■V)fo,=g^ + e^h^, (3.5) 

where 
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The commutator inequality (2.6) and Sobolev embedding theorem imply that 
||.9q|| < C{Me{T)). On the other hand, from the Sobolev embedding theorem and 
the Moscr-typc inequality [23] we get 

I + II /rfvu" • Vii'iii ^ 



|L>"(*'= : Vu"^ 



||/.«!l<i^(||(*^:Vu^)|| 

<C{M,{T)) + C{K{T))\\u%+^. 
Multiplying (3.5) by fa and integrating over [0,t] X with t < T, we obtain 



\\Um' <\\UOW + \\d.u^\\L-iio,ty.w-) / WUrWdT 

\\gc.{T)\\\\fa{r)\\dT + 2e^ f ||/i„(r)||||/„(T)||dT 



<Co + tC{M,{T)) + e^C'iMeiT)), 

where we have used Young's inequality and the embedding H"' ^ L°° for a > d/2. 
The conclusion then follows by adding up these estimates for all jal < s — 1. □ 

The following L^-bound of {q'^, u'^, H*^) can be obtained directly using the energy 
method due to the skew-symmetry of the singular term in the system and the 
special structure of coupling between the magnetic field and fluid velocity. This 
L^-bound is very important in our arguments, since the induction analysis will be 
used to get the desired Sobolev estimates. 

Lemma 3.3. There exist constants Co > and < < and a function C(-) 
independent of e, such that for all t G [0,T], 



(g^ n^Wmf- + ei /* ||Vu^(r)||2rfr < Co + tC{M,{T)). 
Jo 



(3.6) 



Proof. Multiplying (1.16) by q'^, (1.17) by u', and (1.18) by H^, respectively, inte- 
grating over R'', and adding the resulting equations together, we obtain 

{a'dtq\q') + (r^9tu^u^) + {dtW,W) + Ml|Vu^f + (m + A)||divu^|[2 

+ {a'{u' ■ V)q\ q') + {r'{n' ■ V)u% u^) 

= j [(curlH^) X H*^] -vl' dx + j curl (u"^ x H"^) • dx. (3.7) 

Here the singular terms involving 1/e are canceled. Using the identity (1.7) and 
integrating by parts, we immediately obtain that 



j [(curlH^) X W] -u' dx + j curl (u' x H') • W dx 



0. 



In view of the positivity and smoothness of a'^{S'^,eq'^) and r'^ {S'^ , eq'^) , we get 
directly from (1.16), (1.19) and (2.8) that 

\\dtS'\\s-i < C{K{T)), \\edtq^s-i < C{K{T)), (3.8) 

while by the Sobolev embedding theorem, we find that 

\\{dta%dtr')\\L^ < \\{dta\dtr')\l-2 < C{K{T)). (3.9) 

By the definition of Aft (T) and the Sobolev embedding theorem, it is easy to see 
that 

||(Va^Vr^)|U=e <C{M,{T)). 
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Since > 0, 2/i + dA > 0, there exists a positive constant ki such that 
/i||Vu^f + + A)||divu^f > ACillVu^f . 
Thus, from (3.7) wc get that 

{a'q',q') + {r'u\u') + {U%W) + Ki f ||Vu'(r)fdT 

Jo 

+ CiM^T)) f {|<j^(r)p + |u^(r)p + |H^(r)p}dr. (3.10) 
Jo 

Moreover, we have 

Uf + \\u^f < ||(a^)-i|Uoo(aV,9') + ||(r^)-iL=c,(r^u%u^) 

<Co{{a^q\q') + {r'n\u^)), (3.11) 

since a*^ and r"^ are uniformly bounded away from zero. Applying Gronwall's Lemma 
to (3.10), we conclude 

\\{q^,n\W){t)r < CoI|(gg,u5,H^)f exp{iC(A^,(T))}. 

Therefore, the estimate (3.6) follows from an elementary inequality 

e'^'* <l + Ct, 0<t< To, (3.12) 

where Tq is some fixed constant. □ 

Concerning the desired higher order estimates, we cannot directly get them by 
differentiating the system as done in [27], since the coefficients a'^(S"^, eq*^), r'^(S"^, eg*^) 
and b'^{S'^,eq'^) contain two scales S"^ and eq"^. We shall adapt and modify the 
techniques developed in [27] to derive the higher order estimates. Set 

^■'^■-'*)^(*"'"o"'' .-(5.°,-)I„)- ^-(u". 



div 
V 



where 1^ denotes the d x d unit matrix. 

Let LE'{d^) = {E^j-^Lid^) andro(5') =r'^(5'%0). Note that ro (5*=) is smooth, 
positive, and bounded away from zero with respect to each e. First, using Lemma 
3.2 and employing the same analysis as in [27], we have 

Lemma 3.4. There exist constants Ci > 0, K > 0, and a function C{-), depending 
only on Mq, such that for all a £ [1, . . . , s] and t G [0, T], 

IIU^IU < AlL(5,)U^|I,_i + C(|Icurl(roU^)||,_i + ||U'|U_i) (3.13) 

and 

IIU^IU < C{\\{LE^{d.,)Y\J'\\^ + ||curl(roU^)|U-i + ||U^|U-i}, (3.14) 

where C := Ci + tC{M,{T)) + eC{M,{T)). 

We remark that the inequalities (3.13) and (3.14) are similar to the well known 
Helmholtz decomposition, and the estimate on ||S''^(t)||^ in Lemma 3.2 plays a key 
role in the proof of Lemma 3.4. 
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Our next task is to derive a bound on ||{L£;£(9x)}'^U'^||o and ||curl (rou'^)||o--i by 
induction arguments. Let = {i^c (9x)}'^(0, u^)^. We first show the following 
estimate. 

Lemma 3.5. There exist a sujficiently small constant r/i > and two constants 
Co > 0, < ^2 < "^c? a function C(-) from [0,oo) to [0,oo), independent of e, 
such that for all a E [1, . . . , s] and t G [0, T] , 

\\{LMd.)rU'it)r + |£ llVW^f (r)dr 

< Co + tC{K{T)) + ,n /* ||u^(r)||2+idr. (3.15) 



Proof Let IJ% := {LE'{dx)VU'' , a £ {0,...,s}. For simplicity, we set M := 
MeiT), J\f := Me{T), and E {S" , eq") . The case fc = is an immediate 

consequence of Lemma 3.3. It is easy to verify that the operator Lsidx) is bounded 
from H'' to H''-^ for G {1, . . . , s + 1}. Note that the equations (1.16), (1.17) can 
be written as 

{dt + • V)U^ + ^E-^L{dx)V' = E-\3<' + V<=) (3.16) 

with 

(curlH') X y ' div^-" 

For cr > 1; we commute the operator {Le}"^ with (3.16) and multiply the result- 
ing equation by E to infer that 

Eidt + • V)U^ + ii(9x)U^ = E{i„ + + h,), (3.17) 

where 

■.=[dt + u^ ■V,{LeV]U\ 

-MLEViE-'r), 
K -MLEViE-'v'). 

Multiplying (3.17) by and integrating over [0, t] xR'^ with t <T, noticing that 
the singular terms cancel out since L(dx) is skew-adjoint, we use the inequalities 
(3.8) and (3.9), and Cauchy-Schwarz's inequality to deduce that 

{E{t)\J%{t),\J%{t)) <(£;(0)U^(0),U^(0))+C(A^) f\\\5%{T)rdT 

Jo 

+ f ||f,(T)fdT + 2 / / (i?(g, + h,)U^)(r)dr. (3.18) 
Jo Jo Jw 

Following the proof process of Lemma 2.4 in [27], we obtain that 

\\Ut)\\<Cmt))- (3.19) 

Now we estimate the nonlinear term in (3.18) involving g^. We expand g^ as 
follows 

d 

|a|=CT+l 
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d 

i,J = l Ai A2 

:— Bi + B2, 

where 

Ai = {(/3i,..- ,/3fc,7,'5)|l/3i| + --- + |/3fe| + l7l + l'^l <fc + l,0< I7I <fc,|^| <A:}, 

A2 = {l\l = fc + 1 - (lAI + • • • + |/3fc|), (A, • • • 0, 0) e Ai}. 

Since there is no magnetic diffusion in the system, we cannot deal with directly the 
terms involving Bi . Instead, we transform one spatial derivative to . Integrating 
by parts, we have 

d 



d „ 

-EE/ dAE-'Yd",HtH^^\]%dx 

<C{U)+m\\^\r)\\l^i (3.20) 

for sufficiently small constant 771 > 0. 

By virtue of Cauchy-Schwarz's and Sobolev's inequalities, and (3.9), a direct 
computation implies that 



\{EB2V%){T)\dT<C(M). 



(3.21) 



Next, we deal with the term involving the viscosity. Recall that L{dx)^'^ 
(divu', Vg^). Denote 

Li := {a^}-Miv, L2 := {r'Y^W . 
A straightforward computation implies that 



{L1L2} 2 Liu' 



if k is odd; 



if k is even. 



Thus, we induce that 

/ {EKV',){T)dxdT 

An integration by parts gives 
ft 



JR'' 



E{LEV{E-^V')'W'^dxdT. 



f f EL''E{E-^Y')W%dxdT 
Jo Jm^ 

-II /^|VW^|2 + (^ + A)|divW^|2da;dr 
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+ / / E[^E-^l\ + + \)E-^VdiY,{LEY]{0,ll'fV^%dxdT, 

Jo Jw 

where it is easy to verify that 

fc-i 

[E-^A,{LeV] = Y,{LEnE-^A,LE]{LEr-'-\ (3.22) 

i=0 

Noting that LE{dx) = E^^L{dx), we find that 

d 

[E-^A,Le] = -E-^AE-^L{d.,) + J2 B^J^:,,^, 

where Bij = 1, ■ ■ ■ ,d) are the sums of bihnear functions of E~^ and dx{E~^}, 
and Sobolev's inequahties imply that 

||Sz,IU-i<C(AA). 

Thus, we integrate by parts to infer that 

At / / E[E-^A,{LEVmu')'^WldxdT 
Jo Jr'' 

<^ j\\v^%rdr+c{N) j\\wi{T)rdT 

+ f\\mAT)f + \\Hi{T)rdT. 

Jo 

Here i/f is a finite sum of terms of the form 

(9rei)---(a,"'e,)(5M (5>^) 

with \ai\-\ h |a; I + + I7I < cr < s, I7I > 0, and thus |/?| < fc - 1 < s - 1, where 

(ei, . . . , e/), w and denote the coefficients of E~^, Cj and u'^ respectively, with 
Cj taking a form similar to that of Bij . H2 is a finite sum of terms of the form 

with |ai| + ---+|Q;;| + |/3| + |7| < cr + 1 < s + 1, I7I > 1, and thus < cr-l < s-1, 
where {ei, . . . , ei), w and denote the coefficients of E~^, Bij and u*^ respectively. 
Hence, we have 

\\H?r+\\H^r<c{M). 

Similarly, we can show that 



(p + X) f I E[E-^Vdiv,{LEY]{'d,vffV^%dxdT 
Jo Jr'' 



< + ^ 



f \\divW'J^dT + C{N') f \\W'J^ dr + tCiM). 
Jo Jo 



Finally, the above estimates (3.19)-(3.21) and the positivity of E imply (3.15). 

□ 

Next, we derive an estimate for jjcurl (rpu'^) i- Define 

nS\eq^):=l- ''^f\y (3.23) 



INCOMPRESSIBLE LIMIT OF COMPRESSIBLE NON-ISENTROPIC MHD EQUATIONS 13 



Hereafter we denote ro{t) ro{S'^{t)) and /'(i) {S'^ (t) , eq'^ {t)) for notatfonal 
simplicity. 

One can factor out eq'^ in f'^{t). In fact, using Taylor's expansion, one obtains 
that there exists a smooth function g'^{t), such that 

nt) = eg^it) eg^iS^it), eq^{t)), < C(A^,(T)). (3.24) 

Since 

the equations for u'^ are equivalent to 
[dt + {u' ■ V)]{rou') + ivg' =.9'Vg' + (1 - e/)(curlH'^) x H*^ 

+ (1 - 6/)divvI/^ + ^o(g^)" %!(y7^"e) - (3-25) 

We perform the operator curl to the equation (3.25) to obtain that 
[at + (u^-V)](curl(roU^)) 
= [u'^ • V, curl] (rou") + curl [(1 - e/)(curlH') x H'] 
+ curl [(1 - e/)div*'] + [curl,5']Vg' 

Lemma 3.6. There exist constants Cq > 0, < ^3 < /z, a function C(-) from 
[0, 00) to [0, 00) and a sufficiently small constant 772 > 0, such that for all e G (0, 1] 
and all t G [0, T] , 

||{curl(roU^),curlH^}(i)||Li +6 T ||Vcurlu^||2_,dr 

JQ 

<Co+tC{K{T))+n2 I ||u^||'+idT. (3.27) 
Jo 

Proof Set M -.^ K{T) and uj = curl(rou^). Taking d'^ {\a\ < s - 1) on (3.26), 
multiplying the resulting equations by 9"cj, and integrating over [0, t] X with 
t <T, we obtain 

i ||5"w(t) f ^ i ||a"c.(0) f - J\{u^ ■ V)5"a., d^u){r)dr 

+ A[u'-V,a>,a"c.)(r)dr 
Jo 

+ [ (a"{[curl,/]Vg^},5"c.)(r)dr 
Jo 

+ f (a"{[u^- V,curl](roU^)},a"w)(T)dT 
Jo 

+ [ (9"{curl[(l -e/)(curlH') X H']},a"tj)(T)dT 
^0 
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+ / (9"{curl[(l -e/)div^'"]},(9"a;)(T)dr 



1 



(3.28) 



We have to estimate the terms /,;(t) (1 < i < 7) on the right-hand side of (3.28). 
Applying partial integrations, we have 

/i(r) = / |9"c.|2divu' < ||divu^(T)|Uoo|ia"c.(T)f < C(A/-)I|9"c.(r)f , (3.29) 



while for the term I2{t), an application of Cauchy-Schwarz's inequality gives 

|/2(t)| < ||h4r)||, h„(r) ^ ' V,5,>. 

The commutator ha is a sum of terms d^u'^ d'].ijj with multi- indices /3 and 7 satisfying 
+ I7I < s, |/3| > 0, and I7I > 0. Thus, the inequality (2.7) with cr = s - 1 > rf/2 
implies that ||hQ(T)|| < C{N). Hence, we have 



\Ht)\ <C{N) + Wuj{i 



(3.30) 



Noting that ([curl ,g^] a)^ j = aid^^g'^ — ajd^iQ'^ for a = (oi, • • • , a^), the inequality 
(2.7), and the estimate (3.24), we can control the term /3(t) as follows 

l/3(r)| < ||a"{[curl,/]V<7^}|| 

</^|I[eurl,/]V(?^|U_i||5"c^|| 
<A'||V/(T)|U„i||Vg^(r)||,_i||a"c.|| 



< C{N) + \\d"uj{i 



(3.31) 



Similarly, the term /4(t) can be bounded as follows. 

|/4(r)| < A'||9"{K • V,curl](roU<^)}|| W^jW 
</^|I[u^-V,curl](roU^)||,_i||a"w|| 
</^||[u},eurl]a,^.(roU^)|U_i||a"c<.|| 

<C{N) + Wuj{t)\\''. (3.32) 
To bound the term I^ir), we use the Moser-type inequality (see [23]) to deduce 

*^ : Vu^ 



Mr)\ <e'K 



curl 



X u 



<e2x||curlu'||L= 



D 



^fr'^iS")^^ : Vu' 



e2A'||u-^-i(curlu' 



6^(5'% e2q^) 



\d"u;\\ 
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e^K\\D'-^u'\ 



<CW + e'Cmiu^r,^, + ||a"c.(r)f , (3.33) 

where the condition s > 2 + d/2 and the inequahty (2.8) have been used. For the 
term leir), by virtue of (2.1), curl curia = Vdiva — Aa. Thus, we integrate by 
parts to see that 

hir) =(a"{(l -e/)(curlH') x H*^}, ^"{curlcurl (rou^)}) 

= (a" { ( 1 - eg^ ) (curl ) x } , 9" { Vdiv (ro ) - A (ro u*^ ) }) , 

and use Cauchy-Schwarz's inequality and (2.8) to conclude 

|/6(r)| <C(A/-)+(?i||u^(r)||2+i, (3.34) 

where 0i > is a sufficiently small constant independent of e. Next, we deal with 
the term Irir). By the vector identities and integration by parts, we see that there 
exists a sufficiently small 02, such that 

Mr) < - inf{ro(5^)}||Vcurlu^(T)|U_i + C{M) 

+ 02||u^(r)||.+i +eC(A/-)||u^(r)||,+i. (3.35) 

Finally, to estimate ||curlH^||s_i, we apply the operator curl to (1.18) and use 
the vector identity (1.22) to obtain 

at(curlH') • V(curlH') 

= -[curl,u^] • VH^ + curl ((H^ • V)u^ - HMivu^). (3.36) 
By the commutator inequality and Sobolev's inequalities, we find that 

||[curl,u^] • VH^II < C{M) 

and 

||curl((H^ • V)u^ - HMivu^)|| < C(A/') + ^allu^ll'+i 

for sufficiently small constant ^3 > 0. Then, by arguments similar to those used in 
Lemma 3.2, we derive that 

||curlH^||2_^ < Co + <C(A/-) + ^3 f H{T)\\l^,dT. (3.37) 

Jo 

Thus, the lemma follows from adding up the estimates (3.28)-(3.37) for all \a\ < 
s — 1 and choosing constants 9i , 62 and O3 appropriately small. □ 

Next we complete the proof of Theorem 3.1 by the following estimate. 

Lemma 3.7. There exist constants Co > 0, < ^4 < /i, and a function C(-) from 
[0, 00) to [0, 00), such that for all e G (0, 1] and t G [0, T], 

\\{S\q\n\W){t)\\l+U f Wn'W'dr < Co + {t + e)C{M,{T)). (3.38) 

Jo 

Proof. First, from (3.14) we get 

||u^C+i < c{\\Vi{LE^id,)rW% + ||Vcurl(roU^)||^i + ||u^f,}, (3.39) 
where C := Ci+tC{Me{T)) + eC{Me{T)). Moreover, using Lemma 3.2, we obtain 
||u^||^+i < c{\\S7{LE4d.,)rW% + ||V(curlu^)||^_i + ||u^||2}. (3.40) 
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In view of (3.40), there exists a constant K2 such that 

|l|VW^||2+e3||Vcurlu^||2_^ >AC2||u^||^+i-Ci{||V{i^.(a,)rW^||^ 

+ ||V(curlu^)||2_,}-(7||ul2, (3.41) 

where Ci tC{Me{T)) + eC{Me{T)). Now, wc combine the estimates (3.15) and 
(3.27) with (3.41), and use the fact that divH = to conclude that there exists a 
positive constant K3, such that 

\\{LE{d,)rXJ% + \\cm\{rou^)\\l_,+Ks f \\n^\\l+,dx 

Jo 

<Co+tCiMeiT))+C f Wu'WldT 

JQ 

+ C^ 1^ {\\S7{LE^id,)rW^l + \\S7icuAu^)\\l^,]dr 

<Co + {t + e)CiM,iT)) + C f Wu^WldT 

Jo 

for sufficiently small 771 and r/2. Thus by induction, we conclude that 
||{L^,.(5.)rU^||g + ||curl(rou')||^_i 

+ ^3 / HWl+idT <Co + {t + e)C{Me{T)). 

JQ 

Using (3.14) again, we obtain the estimate (3.38) by induction on ct G {0, . . . , s}. □ 

4. Incompressible limit 

In this section, we shall prove the convergence part of Theorem 1.1 by modifying 
the method developed by Metivier and Schochet [27], see also some extensions 
[1,2,26]. 

Proof of the convergence part of Theorem 1.1. The uniform bound (1.29) implies 
that, after extracting a subsequence, one gets the following limits: 

(g^u^H'^) ^ (g,v,H) weakly-* in L°°(0, T; if"(]R'^)). (4.1) 

The equations (1.18) and (1.19) imply that dtS^ andStH^ e C{[Q,T],H^-^{W^)). 
Thus, after further extracting a subsequence, we obtain that, for all s' < s, 

S'^S strongly in C{[0,T], H^Ur")): (4-2) 
H^^H strongly in C([0, T], ifil(M^)), (4.3) 

where the limit H G C([0, T], iJi4(M'^)) n L°°{0,T; Hf^^{R'^)). Similarly, by (3.26) 
and the uniform bound (1.29), we have 

curl {ro{S')u') ^ curl (ro(5)v) strongly in C([0, T], H^l^'iR")) (4.4) 

for aU s' < s, where roiS) = limg^o roiS") := lime_^o r'=(S'% 0). 

In order to obtain the limit system, we need to prove that the convergence in 
(4.1) holds in the strong topology of L'^{0,T; H{^^{R'^)) for aU s' < s. To this end, 
we first show that q = and divv = 0. In fact, from (3.16) we get 

eE%S\eq')dtV' + L{d^)V' -eE%S\ eq')u' ■ VU^ + e(J^ + V^). (4.5) 
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Since 

^;^(5^eg')-^;^(5^0) = O(e), 

we have 

eE'{S', 0)dtU' + L{d.:,)W = eh^ (4.6) 

where is uniformly bounded in C{[0,T], H'-'^iW^)) in view of (1.29). Pass- 
ing to the weak hmit to (4.6), we obtain Vg = and divv — 0. Since q £ 
i°°(0, T; iJ*(M'')), we infer that q — 0. Noticing that the strong compactness 
for the incompressible components by (4.4), it is sufficient to prove the following 
proposition on the acoustic components. 

Proposition 4.1. Suppose that the assumptions in Theorem 1.1 hold, then q'^ 
converges strongly to in L'^{0,T; Hf^^{S.'^)) for all s' < s, and divu' converges 
strongly to in L'^{0,T; H^l'^iR'^)) for all s' < s. 

The proof of Proposition 4.1 is built on the the following dispersive estimates on 
the wave equations obtained by Metivier and Schochct [27] and reformulated in [2]. 

Lemma 4.2 ( [2,27]). LetT > andw" be a bounded sequence inC{[0,T], H'^{W^)), 
such that 

e^dt{a^dtw') -V ■ {b'Vw") ^ c\ 

where c'^ converges to strongly in L^{0, T; L^(R'')). Assume further that, for some 
s > d/2 + 1, the coefficients {a'^,b^) are uniformly bounded in C{[0,T]; H^{R'^)) and 
converges in C{[0,T]; H^^^{M.'^)) to a limit {a,b) satisfying the decay estimate 

\a{x,t)-a\ < Co\x\~^~\ |V,a(a;,t)| < Co\x\~^-\ 

\b{x,t)-b\ < Co\x\-^-', \V^b{x,t)\ < Co\xr^~', 

for some given positive constants a, b, Co and 6. Then the sequence w'^ converges 
toO m i2(o,T;Lf^^(R'^)). 

Proof of Proposition 4-1- We first show that converges strongly to in L^(0, T; 
i?[Qj,(R'')) for all s' < s. An application of the operator e^dt to (1.16) gives 

e^dtia%S\eq')dtq') + edtdivu' = -e^dt{a%S\eq')iu' ■ V)q'}. (4.7) 

Dividing (1.17) by r*^ (S'^ , eq'^) and then taking the operator div to the resulting 
equation, one has 

= -div((u^ • V)u^) + div (-^^^i-^ (curl H^) x H^) 

Subtracting (4.8) from (4.7), we obtain 

e^dt{a%S\eq^)dtqn ~ div[-^^^—^^Vq^) = F%S\ q^ , n\W), (4.9) 

where 

F^(5^g^u^ff) =ediv(— ^l-^(curiff) x h^) 
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+ ediv(-^^^^i-^divv|/(u^)) - ediv((u' • V)u^) 

-e^dt{a\S\eq'){n' ■V)q'}. 

In view of the uniform boundedness of (S*^, q*^, u*^, H*^), the smoothness and pos- 
itivity assumptions on a'^ and r*^, and the convergence of S*', we find that 

F''{S\q\u\W) strongly in L^({),T;L^{M.'^)), 

and the coefficients in (4.9) satisfy the requirements in Lemma 4.2. Therefore, by 
virtue of Lemma 4.2, 

strongly in £^(0, T; if„^(R'')). 

On the other hand, the uniform boundedness of q'^ in C{[Q,T],H'^{W^)) and an 
interpolation argument yield that 

q^^Q strongly in L^{Q,T; H^^^{M!^)) for all s' < s. 

Similarly, we can obtain the convergence of divu*^. □ 

We continue our proof of Theorem 1.1. From Proposition 4.1, we know that 

divu'^divv in L^{Q,T\H^^-'^{W^)). 

Hence, from (4.4) it follows that 

u"^ ^ V in L2(0, T; H^^{W^)) for aU s' < s. 

By (4.2), (4.3) and Proposition 4.2, we obtain 

r^iS^eq^) ^ro{S) in L°°{Q,T;L°°{W^)); 
Vu<^^Vv in L'^{0,T]H^l^^{M.'^))- 
VH^^VH in L'^{Q,T;H^^-^{W^)). 

Passing to the limit in the equations for S*^ and H', we see that the limits S and 
H satisfy 

dtS + (v ■ V)S = 0, 9tH + (v • V)H - (H • V)v = 

in the sense of distributions. Since r'^{S'^, eq'^) — ro{S'^) = 0(e), we have 
ir%S',eq') - roiS^mOtu' + K ' V)u^) ^ 0, 

whence, 

r'iS',eq')idtu' + (u^ • V)u^) = (r'(5^ eg') - ro(5'))(ftu' + (u' • V)u') 

+ 9^(ro(5')u') + (u'-V)(ro(5')u') 

in the sense of distributions. 

Applying the operator curl to the momentum equation (1.17) and taking to the 
limit, we find that 

curl {roiS)idtv + v • Vv) - (curlH) x H - /xAv) = 0. 

Therefore, by the fact that curl V = 0, the limit (S, v, H) satisfies 

r{S, 0)(9tv + (v ■ V)v) - (curlH) x H - /zAv + Vtt = 0, (4.10) 

9tH+(vV)H-(H-V)v = 0, (4.11) 

dtS+{v-\/)S = 0, (4.12) 
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divv = 0, divH = (4.13) 

for some function tt. 

If we employ the same arguments as in the proof of Theorem 1.5 in [27], we 
find that (5,v, H) satisfies the initial conditions (1.31). Moreover, the standard 
iterative method shows that the system (4.10)-(4.13) with initial data (1.31) has a 
unique solution (S**, v*, H*) G C([0, T], iJ*(M'')). Thus, the uniqueness of solutions 
to the limit system (4.10)-(4.13) implies that the above convergence results hold 
for the full sequence of {S'^, 9'^, u*^, H*^). Thus, the proof is completed. □ 



5. Compressible non-isentropic MHD equations with infinite 

Reynolds number 

In the study of magnetohydrodynamics, for some local processes in the cosmic 
system, the effect of the magnetic diffusion will become very important, see [13]. 
Moreover, when the Reynolds number of a fluid is very high and the temperature 
changes very slowly, we can ignore the viscosity and the heat conductivity of the 
fluid in the MHD equations. In such situation, the compressible MHD equations in 



the non-isentropic case take the form: 

+ div(pu) = 0, (5.1) 

dt{pu) + div (pu ® u) + Vp = (curlH) x H, (5.2) 

dt£ + div {u{£' + p)) = div [(u X H) X H + i/H X (curl H)] , (5.3) 

dtH - curl (u X H) = -curl (i^curlH), divH = 0. (5.4) 



As before here p denotes the density, u G M'' (d = 2, 3) the velocity, H e K.'' the 
magnetic field; 8 the total energy given hy £ = £' + |H|^/2 and = p (e + |u|-^/2) 
with e being the internal energy, p|u|^/2 the kinetic energy, and |H|^/2 the magnetic 
energy. The equations of state p ~ p{p, 0) and e = e(p, 0) relate the pressure p and 
the internal energy e to the density p and the temperature 6. The constant v > Q 
is the magnetic diffusivity acting as a magnetic diffusion coefficient of the magnetic 
field. 

Using the Gibbs relation (1.9) and the identities (1.6) and (1.7), the equation of 
energy conservation (5.3) can be replaced by 

dt{pS) + div(pS'u) = ^|curlH|2, (5.5) 
a 

where S denotes the entropy. 

As in Section 1, we reconsider the equations of state as functions of S and p, i.e., 
p — R{S,p) and 9 = Q{S,p) for some positive smooth functions R and 6 defined for 
all S and p > 0, and satisfying dR/dp > 0. Then, by utilizing (1.1) together with 
the constraint divH = 0, the system (5.1), (5.2), (5.4) and (5.5) can be rewritten 
as 

A{S,p){dtP+{u-V)p)+dWu^O, (5.6) 
R{S,p){dtu + (u • V)u) + Vp = (curlH) x H, (5.7) 
^tH - curl (u X H) = -curl (z/ curl H), divH = 0, (5.8) 
R{S,p)Q{S,p){dtS + (u • V)5) = i/|curlH|2, (5.9) 
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where A{S,p) is defined by (1.15). By introducing the dimensionless parameter e, 
and making the following changes of variables: 

p(x, t) = p'(.T, et), S{x, t) = S'ix, et), 

u(a-, t) = eu^(a;, et), 'ii{x,t) — tiiV{x,et), v = t , 

and p'^{x,et) = pe'^'^ (.^'<^t) for some positive constant p, the system (5.6)-(5.9) can 
be rewritten as 

a'(5^ eq'){dtq'' + {u' ■ V)<?') + ^divu' = 0, (5.10) 

e 

r'{S',eq')idtu' + (u^ • V)u^) + ^Vg^ = (curlff ) x H^ (5.11) 

e 

dtH' - curl (u' X H') - ^"AH' = 0, divH"^ = 0, (5.12) 

b'{S',eq'){dtS' + (u' ■ \/)S') = eVicurlH'p, (5.13) 

where we have used the identity curl curl H'^ = VdivH'^ — AH*^, the constraint 
divH^ = 0, and the abbreviations (1.20) and (1.21). 

Formally, we obtain from (5.10) and (5.11) that S/q'^ and divu^ = as 
e — > 0. Applying the operator curl to (5.11), using the fact that curlV = 0, and 
letting e — 0, we obtain 

curl {r{S, 0){^t^r + v Vv) - (curlH) x H) = 0, 

where we have assumed that (S"^, g*^, u'^, H^) and r'^(S"^, eg*^) converge to {S, 0, v, H) 
and r(S, 0) in some sense, respectively. Finally, Letting /i*^ — > > and applying 
the identity (1.22), we expect to get the following incompressible MHD equations: 

r(5,0)(atv+(v-V)v)-(curlH) x H + Vtt 0, (5.14) 

9tH+(vV)H-(H-V)v-^AH = 0, (5.15) 

dtS+{vW)S = 0, (5.16) 

divv = 0, divH = (5.17) 

for some function tt. 

We supplement the system (1.16)-(1.19) with initial conditions 

(5^g^u^H^)|t^o = (5^,g^,u^,H^). (5.18) 
The main result of this section reads as follows. 

Theorem 5.1. Let s > d/2 + 2 be an integer and — /i > 0. For any constant 
Mq > 0, there is a positive constant T = T[Mq), such that for all e G (0,1] and 
any initial data (^g, qg, u^, Hg) satisfying 

||(5^,qS,u^,H^)||a=(K.) < A./o, (5.19) 

the Cauchy problem (5.10)-(5.13), (5.18) has a unique solution (S**^, q'^, u*^, H') £ 
C0([0,r],7?''(R'')), satisfying that for all e G (0,1] and t e [0,T], 

||(S",g%u%H")(i)||jj=(K'') < ^ for some constant N ^ N{Mo) > 0. (5.20) 

Moreover, suppose further that the initial data (5q, gg, u^, Hq) converge strongly in 
if^(IR'') to (5*0, 0, vg. Ho) and 5q decays sufficiently rapidly at infinity in the sense 
that 

\S'oix)'S\<No\x\-'-\ \yS'oi^)\<No\x\-^-' (5.21) 
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for all e £ (0,1] and some positive constants S_, Nq and l. Then, (5*^, g*^, u*^, H*^) 
converges_ weakly in L°°{Q,T;H''{W^)) and strongly in L'^{Q,T;Hf^^{W^)) to a limit 
(5, 0,v,H) for all s' < s. Moreover, (S'jVjH) is the unique solution in C([0,r], 
H-''{R'^)) to the system (5.14) -(5.17) with initial data (5,v,H)|(=o = (^o, Wq, Hq), 
where wq G H^{M.'^) is determined by 

divwo = 0, curl (r(S'o, 0)wo) = curl (r(S'o,0)vo), r(5o, 0) lim r"(5^, 0). (5.22) 

The function tt G C([0,T] x R'^) satisfies Vtt G C{[0,T], H'-\R'^)). 

Sketch of the proof of Theorem 5.1. As explained before, the main step is to estab- 
lish the uniform estimate (5.20). For this purpose, we define A4e{T) as follows 

M,{T) -.^KiTf + r ||H^||^+idT, (5.23) 
Jo 

where Afe{T) is defined by (3.3). By arguments similar to those used in the proof 
of Theorem 3.1, one can obtain the desired estimate. Indeed, the arguments are 
easier since one can use the magnetic diffusion term to control the terms involving 
H in the momentum equations, and therefore we omit the details here. □ 
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